Histone Modifications at Human Enhancers Reflect Global Cell-Type-Specific Gene Expression

The human body is composed of diverse cell types with distinct functions. Although it is known that lineage specification depends on cell-specific gene expression, which in turn is driven by promoters, enhancers, insulators and other cis-regulatory DNA sequences for each gene1, 2, 3, the relative roles of these regulatory elements in this process are not clear. We have previously developed a chromatin-immunoprecipitation-based microarray method (ChIP-chip) to locate promoters, enhancers and insulators in the human genome4, 5, 6. Here we use the same approach to identify these elements in multiple cell types and investigate their roles in cell-type-specific gene expression. We observed that the chromatin state at promoters and CTCF-binding at insulators is largely invariant across diverse cell types. In contrast, enhancers are marked with highly cell-type-specific histone modification patterns, strongly correlate to cell-type-specific gene expression programs on a global scale, and are functionally active in a cell-type-specific manner. Our results define over 55,000 potential transcriptional enhancers in the human genome, significantly expanding the current catalogue of human enhancers and highlighting the role of these elements in cell-type-specific gene expression

Discovery and characterization of chromatin states for systematic annotation of the human genome

A plethora of epigenetic modifications have been described in the human genome and shown to play diverse roles in gene regulation, cellular differentiation and the onset of disease. Although individual modifications have been linked to the activity levels of various genetic functional elements, their combinatorial patterns are still unresolved and their potential for systematic de novo genome annotation remains untapped. Here, we use a multivariate Hidden Markov Model to reveal 'chromatin states' in human T cells, based on recurrent and spatially coherent combinations of chromatin marks. We define 51 distinct chromatin states, including promoter-associated, transcription-associated, active intergenic, large-scale repressed and repeat-associated states. Each chromatin state shows specific enrichments in functional annotations, sequence motifs and specific experimentally observed characteristics, suggesting distinct biological roles. This approach provides a complementary functional annotation of the human genome that reveals the genome-wide locations of diverse classes of epigenetic function.National Science Foundation (U.S.). (Award 0905968)National Human Genome Research Institute (U.S.) (Award U54-HG004570)National Human Genome Research Institute (U.S.) (Award RC1-HG005334

Chromatin States Accurately Classify Cell Differentiation Stages

Gene expression is controlled by the concerted interactions between transcription factors and chromatin regulators. While recent studies have identified global chromatin state changes across cell-types, it remains unclear to what extent these changes are co-regulated during cell-differentiation. Here we present a comprehensive computational analysis by assembling a large dataset containing genome-wide occupancy information of 5 histone modifications in 27 human cell lines (including 24 normal and 3 cancer cell lines) obtained from the public domain, followed by independent analysis at three different representations. We classified the differentiation stage of a cell-type based on its genome-wide pattern of chromatin states, and found that our method was able to identify normal cell lines with nearly 100% accuracy. We then applied our model to classify the cancer cell lines and found that each can be unequivocally classified as differentiated cells. The differences can be in part explained by the differential activities of three regulatory modules associated with embryonic stem cells. We also found that the “hotspot” genes, whose chromatin states change dynamically in accordance to the differentiation stage, are not randomly distributed across the genome but tend to be embedded in multi-gene chromatin domains, and that specialized gene clusters tend to be embedded in stably occupied domains

Global analysis of DNA methylation in early-stage liver fibrosis

<p>Abstract</p> <p>Background</p> <p>Liver fibrosis is caused by chemicals or viral infection. The progression of liver fibrosis results in hepatocellular carcinogenesis in later stages. Recent studies have revealed the importance of DNA hypermethylation in the progression of liver fibrosis to hepatocellular carcinoma (HCC). However, the importance of DNA methylation in the early-stage liver fibrosis remains unclear.</p> <p>Methods</p> <p>To address this issue, we used a pathological mouse model of early-stage liver fibrosis that was induced by treatment with carbon tetrachloride (CCl₄) for 2 weeks and performed a genome-wide analysis of DNA methylation status. This global analysis of DNA methylation was performed using a combination of methyl-binding protein (MBP)-based high throughput sequencing (MBP-seq) and bioinformatic tools, IPA and Oncomine. To confirm functional aspect of MBP-seq data, we complementary used biochemical methods, such as bisulfite modification and <it>in-vitro</it>-methylation assays.</p> <p>Results</p> <p>The genome-wide analysis revealed that DNA methylation status was reduced throughout the genome because of CCl₄treatment in the early-stage liver fibrosis. Bioinformatic and biochemical analyses revealed that a gene associated with fibrosis, <it>secreted phosphoprotein 1 </it>(<it>Spp1</it>), which induces inflammation, was hypomethylated and its expression was up-regulated. These results suggest that DNA hypomethylation of the genes responsible for fibrosis may precede the onset of liver fibrosis. Moreover, <it>Spp1 </it>is also known to enhance tumor development. Using the web-based database, we revealed that <it>Spp1 </it>expression is increased in HCC.</p> <p>Conclusions</p> <p>Our study suggests that hypomethylation is crucial for the onset of and in the progression of liver fibrosis to HCC. The elucidation of this change in methylation status from the onset of fibrosis and subsequent progression to HCC may lead to a new clinical diagnosis.</p

DNA fragments binding CTCF in vitro and in vivo are capable of blocking enhancer activity

<p>Abstract</p> <p>Background</p> <p>Earlier we identified ten 100-300-bp long CTCF-binding DNA fragments selected earlier from a 1-Mb human chromosome 19 region. Here the positive-negative selection technique was used to check the ability of CTCF-binding human genomic fragments to block enhancer-promoter interaction when inserted into the genome.</p> <p>Results</p> <p>Ten CTCF-binding DNA fragments were inserted between the CMV enhancer and CMV minimal promoter driving the herpes simplex virus thymidine kinase (HSV<it>-tk</it>) gene in a vector expressing also the <it>neo</it>^Rgene under a separate promoter. The constructs were then integrated into the genome of CHO cells, and the cells resistant to neomycin and ganciclovir (positive-negative selection) were picked up, and their DNAs were PCR analyzed to confirm the presence of the fragments between the enhancer and promoter in both orientations.</p> <p>Conclusions</p> <p>We demonstrated that all sequences identified by their CTCF binding both <it>in vitro </it>and <it>in vivo </it>had enhancer-blocking activity when inserted between the CMV minimal promoter and enhancer in stably transfected CHO cells.</p

Complex exon-intron marking by histone modifications is not determined solely by nucleosome distribution

It has recently been shown that nucleosome distribution, histone modifications and RNA polymerase II (Pol II) occupancy show preferential association with exons (“exon-intron marking”), linking chromatin structure and function to co-transcriptional splicing in a variety of eukaryotes. Previous ChIP-sequencing studies suggested that these marking patterns reflect the nucleosomal landscape. By analyzing ChIP-chip datasets across the human genome in three cell types, we have found that this marking system is far more complex than previously observed. We show here that a range of histone modifications and Pol II are preferentially associated with exons. However, there is noticeable cell-type specificity in the degree of exon marking by histone modifications and, surprisingly, this is also reflected in some histone modifications patterns showing biases towards introns. Exon-intron marking is laid down in the absence of transcription on silent genes, with some marking biases changing or becoming reversed for genes expressed at different levels. Furthermore, the relationship of this marking system with splicing is not simple, with only some histone modifications reflecting exon usage/inclusion, while others mirror patterns of exon exclusion. By examining nucleosomal distributions in all three cell types, we demonstrate that these histone modification patterns cannot solely be accounted for by differences in nucleosome levels between exons and introns. In addition, because of inherent differences between ChIP-chip array and ChIP-sequencing approaches, these platforms report different nucleosome distribution patterns across the human genome. Our findings confound existing views and point to active cellular mechanisms which dynamically regulate histone modification levels and account for exon-intron marking. We believe that these histone modification patterns provide links between chromatin accessibility, Pol II movement and co-transcriptional splicing

Discovery and Annotation of Functional Chromatin Signatures in the Human Genome

Transcriptional regulation in human cells is a complex process involving a multitude of regulatory elements encoded by the genome. Recent studies have shown that distinct chromatin signatures mark a variety of functional genomic elements and that subtle variations of these signatures mark elements with different functions. To identify novel chromatin signatures in the human genome, we apply a de novo pattern-finding algorithm to genome-wide maps of histone modifications. We recover previously known chromatin signatures associated with promoters and enhancers. We also observe several chromatin signatures with strong enrichment of H3K36me3 marking exons. Closer examination reveals that H3K36me3 is found on well-positioned nucleosomes at exon 5′ ends, and that this modification is a global mark of exon expression that also correlates with alternative splicing. Additionally, we observe strong enrichment of H2BK5me1 and H4K20me1 at highly expressed exons near the 5′ end, in contrast to the opposite distribution of H3K36me3-marked exons. Finally, we also recover frequently occurring chromatin signatures displaying enrichment of repressive histone modifications. These signatures mark distinct repeat sequences and are associated with distinct modes of gene repression. Together, these results highlight the rich information embedded in the human epigenome and underscore its value in studying gene regulation

Predicting cell types and genetic variations contributing to disease by combining GWAS and epigenetic data

Genome-wide association studies (GWASs) identify single nucleotide polymorphisms (SNPs) that are enriched in individuals suffering from a given disease. Most disease-associated SNPs fall into non-coding regions, so that it is not straightforward to infer phenotype or function; moreover, many SNPs are in tight genetic linkage, so that a SNP identified as associated with a particular disease may not itself be causal, but rather signify the presence of a linked SNP that is functionally relevant to disease pathogenesis. Here, we present an analysis method that takes advantage of the recent rapid accumulation of epigenomics data to address these problems for some SNPs. Using asthma as a prototypic example; we show that non-coding disease-associated SNPs are enriched in genomic regions that function as regulators of transcription, such as enhancers and promoters. Identifying enhancers based on the presence of the histone modification marks such as H3K4me1 in different cell types, we show that the location of enhancers is highly cell-type specific. We use these findings to predict which SNPs are likely to be directly contributing to disease based on their presence in regulatory regions, and in which cell types their effect is expected to be detectable. Moreover, we can also predict which cell types contribute to a disease based on overlap of the disease-associated SNPs with the locations of enhancers present in a given cell type. Finally, we suggest that it will be possible to re-analyze GWAS studies with much higher power by limiting the SNPs considered to those in coding or regulatory regions of cell types relevant to a given disease

Relationship between Gene Body DNA Methylation and Intragenic H3K9me3 and H3K36me3 Chromatin Marks

To elucidate the relationship between intragenic DNA methylation and chromatin marks, we performed epigenetic profiling of chromosome 19 in human bronchial epithelial cells (HBEC) and in the colorectal cancer cell line HCT116 as well as its counterpart with double knockout of DNMT1 and DNMT3B (HCT116-DKO). Analysis of H3K36me3 profiles indicated that this intragenic mark of active genes is associated with two categories of genes: (i) genes with low CpG density and H3K9me3 in the gene body or (ii) genes with high CpG density and DNA methylation in the gene body. We observed that a combination of low CpG density in gene bodies together with H3K9me3 and H3K36me3 occupancy is a specific epigenetic feature of zinc finger (ZNF) genes, which comprise 90% of all genes carrying both histone marks on chromosome 19. For genes with high intragenic CpG density, transcription and H3K36me3 occupancy were not changed in conditions of partial or intensive loss of DNA methylation in gene bodies. siRNA knockdown of SETD2, the major histone methyltransferase responsible for production of H3K36me3, did not reduce DNA methylation in gene bodies. Our study suggests that the H3K36me3 and DNA methylation marks in gene bodies are established largely independently of each other and points to similar functional roles of intragenic DNA methylation and intragenic H3K9me3 for CpG-rich and CpG-poor genes, respectively

Enrichment analysis of Alu elements with different spatial chromatin proximity in the human genome

Transposable elements (TEs) have no longer been totally considered as “junk DNA” for quite a time since the continual discoveries of their multifunctional roles in eukaryote genomes. As one of the most important and abundant TEs that still active in human genome, Alu, a SINE family, has demonstrated its indispensable regulatory functions at sequence level, but its spatial roles are still unclear. Technologies based on 3C(chromosomeconformation capture) have revealed the mysterious three-dimensional structure of chromatin, and make it possible to study the distal chromatin interaction in the genome. To find the role TE playing in distal regulation in human genome, we compiled the new released Hi-C data, TE annotation, histone marker annotations, and the genome-wide methylation data to operate correlation analysis, and found that the density of Alu elements showed a strong positive correlation with the level of chromatin interactions (hESC: r=0.9, P<2.2×1016; IMR90 fibroblasts: r = 0.94, P < 2.2 × 1016) and also have a significant positive correlation withsomeremote functional DNA elements like enhancers and promoters (Enhancer: hESC: r=0.997, P=2.3×10−4; IMR90: r=0.934, P=2×10−2; Promoter: hESC: r = 0.995, P = 3.8 × 10−4; IMR90: r = 0.996, P = 3.2 × 10−4). Further investigation involving GC content and methylation status showed the GC content of Alu covered sequences shared a similar pattern with that of the overall sequence, suggesting that Alu elements also function as the GC nucleotide and CpG site provider. In all, our results suggest that the Alu elements may act as an alternative parameter to evaluate the Hi-C data, which is confirmed by the correlation analysis of Alu elements and histone markers. Moreover, the GC-rich Alu sequence can bring high GC content and methylation flexibility to the regions with more distal chromatin contact, regulating the transcription of tissue-specific genes